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Abstract 

In the years from 2001 to 2006, Guseinov and his coworkers pubhshed 40 arti- 
cles on the derivation and appUcation of one-range addition theorems. In E. J. 
Weniger, Extended Comment on "One-Range Addition Theorems for Coulomb 
Interaction Potential and Its Derivatives" by I. I. Guseinov (Chem. Phys. Vol. 
309 (2005), pp. 209 - 213), arXiv:0704.1088v2 [math-ph], it was argued that 
Guseinov's treatment of one-range addition theorems is at best questionable 
and in some cases fundamentally flawed. In I. I. Guseinov, Extended Rejoin- 
der to "Extended Comment on "One- Range Addition Theorems for Coulomb 
Interaction Potential and Its Derivatives" by I. I. Guseinov (Chem. Phys. and 
Vol. 309 (2005)", pp. 209-213), arXiv:0706.0975v2 [physics.chem-ph], these 
claims were disputed. To clarify the situation, the most serious mathematical 
flaws in Guseinov's treatment of one-range addition theorems are discussed in 
more depth. 
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1 Introduction 

The efficient and reliable evaluation of multicenter integrals is among the oldest 
mathematical and computational problems of molecular electronic structure theory. 
In spite of heroic efforts, the mathematical and computational problems, that occur 
in this context, are not yet solved in a completely satisfactory way, and there is still 
a considerable amount of research going on. Particularly difficult are multicenter 
integrals of the physically better motivated exponentially decaying functions, whose 
efficient and reliable evaluation is - in spite of all the mathematical and computational 
advances of recent years - still very difficult. 

From a methodological point of view, research on multicenter integrals is essen- 
tially mathematical in nature, although relatively few mathematicians have been 
involved in this research. In my opinion, this is quite deplorable. But within math- 
ematics, research on multicenter integrals is highly interdisciplinary. Absolutely es- 
sential is a good knowledge of special function theory and of classical analysis: The 
derivation of explicit expressions for multicenter integrals is to a large extent some 
kind of 19th century mathematics. However, in order to succeed we also need modern 
mathematical concepts as for example Hilbert spaces, approximation theory, gener- 
alized functions, and angular momentum theory. 

It is the ultimate goal of research on multicenter integrals to produce computer 
code that permits an efficient and reliable evaluation of these integrals. Accordingly, 
a good knowledge of sophisticated numerical techniques is absolutely indispensable. 

Research on multicenter integrals is difficult, and there are many chances of mak- 
ing errors. Firstly, there are errors that violate basic mathematical principles, which 
could be called first- order errors. Secondly, there are mathematically correct manip- 
ulations and/or deductions, which lead to computer code that is either hopelessly 
inefficient or unreliable and which could be called second- order errors. 

Multicenter integrals are difficult to evaluate because the integration variables 
occur in unseparated form. Principal mathematical tools, that can accomplish such 
a separation of variables, are so-called addition theorems. These are expansions of 
a given function f{r ± r') with r,r' G M.^ in products of other functions that only 
depend on either r or r'. 

In the articles [Il[2l[3lia[5l[6l[7l[8l[9l[T0l[IIl[l2l[l3l[Tl[T^ 
ini[22l[23l[2a[25l[26l[27l[28l[29l[3a Guseinov 
and coworkers derived so-called one-range addition theorems for a variety of different 
functions and applied them for the evaluation of multicenter integrals. 

In my Comment [H], I presented a very detailed criticism of Guseinov's work 
on one-range addition theorems and showed that both ffist- and second-order errors 
occur quite abundantly. I also criticized in [41j that Guseinov's articles are often re- 
markably similar and that they do not always give due credit to the previous work of 
others. Moreover, I made several suggestions how Guseinov's flawed approach could 
be improved: In [HI Section 7] I suggested to employ advanced mathematical con- 
cepts from the theory of generalized functions in order to give a meaning to divergent 
one-range addition theorems, or to use nonlinear sequence transformations (see for 
example [I2ll33] and references therein), either to sum divergent series, whose occur- 
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rence was apparently either overlooked or ignored, or to accelerate the convergence 
of series expansions for multicenter integrals, whose convergence need not be rapid. 

Apparently, my Comment [41j did not impress Guseinov too much, who recently 
wrote a Rejoinder [44J to my Comment [H] In his Abstract, Guseinov stated 

The concrete criticism raised in Weniger 's comment against our papers 
actually touches a very minor aspect of the works that are not relevant 
at all for the conclusions, which are made. 

and 

All claims of inconsistencies and flaws in the theoretical framework are 
rejected as unfounded. This rejoinder paper contains all of the answers 
to Weniger 's comments. 

This is wrong. In pH], Guseinov did not address at all the most serious and conse- 
quential first-order errors, which I had criticized in my Comment [H] and which raise 
serious doubts on the mathematical soundness of Guseinov's treatment of one-range 
addition theorems. Guseinov's questionable attitude towards mathematical rigor be- 
comes also evident in his even more recent preprints P5| B6]. in which he proceeds in 
the same style and spirit as in his previous articles on one-range addition theorems 
mentioned above, completely ignoring my criticism. 

In this Replay to Guseinov's Rejoinder |44j, I discuss once more and in more 
depth the most important mathematical flaws of Guseinov's work. In contrast to my 
earlier and longer Comment [H], I concentrate entirely on first-order errors. 

In [U E] , Guseinov had derived one-range addition theorems for Slater-type func- 
tions 

x'JMr) = iPrf-'^-'e-f'^y^iPr), /5 > , (1.1) 

with in general nonintegral principal quantum numbers G M \ N by expanding 
them in terms his functions r) defined by (ID. 10 . which are complete and 

orthonormal in the weighted Hilbert spaces L^^. (M^) defined by (ID. 51) (for details, see 
[m Sections 4 - 6]). Here, 3^1^ (/3r) is a regular solid harmonic (compare [H], Eq. 
(A.2)])). 

As long as the principal quantum numbers A^ are not too negative. Slater-type 
functions Xnl belong for k = —1, 0, 1, 2, . . . to the weighted Hilbert spaces L^fc(M^) 
defined by fID.Sp . which was implicitly used by Guseinov. In this case, the one-range 



addition theorems derived in |ll [3] exist and converge in the mean with respect to 
the norms flD.4p of these Hilbert spaces. 

Guseinov's addition theorems for Slater-type functions Xnl yield for N = L = 
M = the corresponding addition theorems for the Yukawa potential exp(— /?r)/r. 
Guseinov [17] derived his one-range addition theorems for the Coulomb potential 
by exploiting the obvious relationship 1/r = \imj3^Qexp{—f3r)/r in his one-range 
addition theorems for the Yukawa potential. 

At first sight, it may look like a good idea to derive one-range addition theorems 
for the Coulomb potential by performing the comparatively simple limit /3 — > in 
the one-range addition theorems for the Yukawa potential. This pragmatic approach 
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seems to permit an efficient utilization of the available information and of previous 
work. Unfortunately, the situation is much more complicated. 

As shown in Section [21 the limit /5 — >■ in mathematically well defined integrals 
containing the Yukawa potential is not necessarily continuous and does not always 
lead to a finite result. So, whenever we try to find an expression for an integral con- 
taining the Coulomb potential as the limiting case of the expression for an analogous 
integral containing the Yukawa potential, we have to be cautious and be prepared 
for complications. As discussed in Sections [3] and HI similar complications can occur 
also in orthogonal expansions of the Yukawa potential. 

A very serious weakness of Guseinov's earlier work on one-range addition theorems 
is that he completely ignored the obvious fact that orthogonal expansions can diverge, 
and he still does this in his recent Rejoinder [H] and in his even more recent preprints 
SS]. On p. 3 of his Rejoinder [H], Guseinov states: 

The essential facts of Hilbert space and approximation theory as weU as 
all questions of convergence and existence have been taken into account by 
Guseinov and his coworkers in the context of one-range addition theorems 
and multicenter integrals. Thus, the Guseinov's treatment of one-range 
addition theorems is not questionable, and is fundamentally flawless from 
a mathematical point of view. 

Again, this is wrong. Let us assume that H is a Hilbert space, and that {v^nj^o ^ 
function set that is complete and orthonormal in ?-^. If / G 7i, then standard Hilbert 
space theory tells us that the expansion ^^o(v^"l/)'^" converges to / in the mean 
with respect to the norm of Ti. If, however, / ^ 7i, then it follows - as shown in 
Section [3]- from the Riesz-Fischer Theorem (see for example Theorem 7.43 on p. 
191]) that the formal orthogonal expansion / = ^^o(v'"l/)'^" diverges in the mean 
with respect to the norm of 7i. 

Since Guseinov ignores the fact that orthogonal expansions can diverge, he fails 
to take into account that divergent orthogonal expansions cannot be treated like 
convergent orthogonal expansions (obviously, the divergence of an orthogonal series 
can do a lot of harm in integrals). So, Guseinov never wonders whether and under 
which conditions divergent expansions can safely be used in multicenter integrals, 
although this is by no means obvious. 

As shown in Section [31 it is nevertheless possible to use divergent orthogonal 
expansions in inner products or other functionals - in our case usually multicenter 
integrals - in a mathematically meaningful way. The key is that the functionals, in 
which these expansions are to be used, have to satisfy additional and possibly very 
restrictive regularity conditions. This follows from the Riesz Representation Theorem 
(see for example |47|, Theorem 7.60 on p. 199]). Accordingly, divergent orthogonal 
expansions are essentially generalized functions in the sense of Schwartz [48j that 
can - in spite of their divergence - converge weakly when used in suitably restricted 
functionals. 

When Guseinov derived in [T7] one-range addition theorems for the Coulomb 
potential by considering the limit — *• in one-range addition theorems for the 
Yukawa potential exp(— /3r)/r, he overlooked some very consequential facts. 
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The Yukawa potential belongs for = 0, 1, 2, . . . , to the weighted Hilbert space 
L^fc(M^)l defined by ( ID. 511 . but not for k = —1 (this obvious fact was apparently 
overlooked by Guseinov). Accordingly, Guseinov's one-range addition theorems for 
the Yukawa potential converge for = 0, 1, 2, . . . in the mean with respect to the norm 
( ID.4I) of L^fe(M^), but not for k = —1 (see also Appendix [El where the convergence of 
the one-center limit of these addition theorems is analyzed). 

As discussed in Section Hj the Coulomb potential does not belong to any of the 
Hilbert spaces L^fe(M^) which Guseinov had implicitly used. This is quite consequen- 
tial: Guseinov's one-range addition theorems for the Coulomb potential diverge for all 
k = —1, 0, 1, 2, . . . in the mean with respect to the norm (10.40 of L^j.(M'^), although 
they were derived by a limiting procedure from the corresponding addition theorems 
of the Yukawa potential, which converge at least for /c = 0, 1, 2, . . . in the mean. This 
is another example that the limit /? — *• in expressions involving the Yukawa poten- 
tial exp(— /5r)/r need not be continuous and does not necessarily produce something 
finite. 

Section HI shows conclusively that it is impossible to construct one-range addition 
theorems for the Coulomb potential that converge in the mean with respect to the 
norms of the weighted Hilbert spaces L^fc(M^). More advanced mathematical concepts 
such as the theory of generalized functions or possibly also powerful numerical tech- 
niques for the summation of divergent series (see for example [l2l B3] and references 
therein) are needed to give divergent orthogonal expansions of that kind any meaning 
beyond purely formal expansions. None of these things are discussed in Guseinov's ar- 
ticles on one-range addition theorems PElEliaElEllTllHllHliTnillllliaiTHliTliTSl 
ini [181 [201 lai E2l [231 [22 [251 [261 [23 
or in his recent Rejoinder [H]. 

Unfortunately, this is not yet the end of Guseinov's grave first-order errors in 
his treatment of one-range addition theorems. As discussed in Section [3, Guseinov 
preferred to replace the complete and orthonormal expansion functions ^\l/™^(/5, r) of 
his one-range addition theorems for Slater-type functions by nonorthogonal Slater- 
type functions with integral principal quantum numbers and to rearrange the order of 
summations of the resulting expansions. In this way, Guseinov obtained expansions 
of Slater-type functions lil^^ r ±r') with in general nonintegral principal quantum 
numbers N e R\N in terms of Slater-type functions X^iil^-, r) with integral principal 
quantum numbers n G N located at a different center. 

As is well known. Slater-type functions are complete in the Hilbert spaces im- 
plicitly used by Guseinov, but not orthogonal. This is very consequential. It is 
extensively documented both in the mathematical literature (see for example [IHl 
Theorem 10 on p. 54] or [SUl Section 1.4]) as well as in the literature on electronic 
structure calculations [511 [521 [531 [SH [551 [56] ) that the existence of expansions in terms 
of nonorthogonal function sets is not guaranteed in the case of essentially arbitrary 
functions. Such an expansions may or may not exist. Thus, as already emphasized 
in [m Section 6]), Guseinov's approach is dangerous and potentially disastrous and 
the validity of his rearranged addition theorems has to be checked explicitly. 

As discussed in Section [5l Guseinov disagreed in his Rejoinder [HI p. 7] with my 
conclusions, and claimed instead that the validity of his approach follows from Eq. 
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(3.11) of his Rejoinder. It is easy to sliow that Guseinov's reasoning is superficial 
and that it is indeed necessary to analyze whether Guseinov's rearranged addition 
theorems exist or not. 

There is the practical problem that one-range additions theorems for exponentially 
decaying functions are fairly complicated mathematical objects. Accordingly, explicit 
proofs of their convergence or divergence are very difficult and would most likely 
require a considerable amount of time and effort. Fortunately, at least some insight 
can be gained by analyzing instead the much simpler one-center limits of Guseinov's 
one-range addition theorems, although this approach does not answer all questions 
of interest. 

In Section O and also in [HI Section 6], it is shown that it is impossible to rear- 
range the one-center limits of Guseinov's one-range addition theorems for Slater-type 
functions if the principal quantum number is nonintegral, iV G M \ N. Guseinov's 
rearrangements do not lead to divergent series in the usual sense, but to power series 
with series coefficients that are for all but a finite number of indices infinite. While 
I can probably claim with some confidence that I have a lot of experience with the 
summation of divergent series (see for example |l2l |33[ [57] and references therein), 

1 nevertheless must admit that I have not the slightest idea what to do with power 
series with an infinite number of infinite terms. 

This article is concluded by a Summary in Section ([6|). For the convenience of 
the readers, the most important conventions and definitions of this Reply are listed 
in Appendices ([A]) - (|D]). Finally, there is Appendix ([E|) analyzing the convergence 
of the one-center limit of Guseinov's one-range addition theorems for the Yukawa 
potential in the weighted Hilbert space L^fe(M^). 

2 On the Continuity of Limits in Integrals 

Because of its exponential decay, the Yukawa potential exp(— /3r)/r is in many re- 
spects a much more convenient mathematical object than the closely related Coulomb 
potential 1/r. This is particularly true for integrals over the whole three-dimensional 
space as they occur in atomic or molecular electronic structure calculations. 

Consequently, it is an obvious idea to derive explicit expressions for integrals 
involving the Coulomb potential by performing the limit /3 — *• in explicit expressions 
for analogous integrals involving the more convenient Yukawa potential. 

Often, this indirect approach is very effective. However, it is no panacea. More- 
over, it can easily lead to problems: The limiting process 1/r = lim^^o exp(— /5r)/r 
is not necessarily continuous in integrals, and it is not guaranteed that it produces a 
finite result. 

These possible problems can be illuminated easily by considering integrals of the 
following kind: 

J{f;P) = I [f{r)Y'-^^tMd^r (2.1) 
As usual, integration extends over the whole R^. 
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For P > 0, the Yukawa potential belongs to the Hilbert space L^(]R^) of square 
integrable functions defined by ( IB. 51) . If we also have / G L^(M^), the integral (12. ip 
is a special case of the inner products (]B.3p . and it is finite. This follows at once form 
the Cauchy-Schwarz Inequality (see for example [17|, Theorem 7.7 on p. 177]) which 
can be expressed as follows: 

\{f\9)\' < {f\f){9\9) = ll/flkf ■ (2.2) 

If we perform the limit /3 in (12. ip . we formally obtain integrals of the following 
kind: 

/C(/) = I [/(r)]*idV (2.3) 

Since the Coulomb potential does not belong to L^(R'^) or to any of those Hilbert 
spaces, which are considered in this Reply and which all involve an integration over 
the whole M.^, the Cauchy-Schwarz Inequality cannot be used to guarantee the exis- 
tence of /C(/) for arbitrary / G L^(M^). Instead, the integral )C{f) makes sense only 
if / belongs to a suitably restricted (proper) subset of L^(]R^) pr to other function 
spaces. 

The possible discontinuity of the limit j3 ^ also becomes evident in the six- 
dimensional integrals 

C{f,g) = j j [f{r)Yy^^9{r')drdr' (2.4) 

and 

yU.9;P) = [/(r')]* '^P^^/_|'^^~^^'^ (rOdrdr-, (2.5) 

which describe the interaction of two charge densities /, (7 : M'^ — > C via the Coulomb 
and the Yukawa potential, respectively, and which involve an integration over the 
whole six-dimensional space x M'^. 

Obviously, we have C{f,g) = lim^^o 3^(/) fi'S /?), but it would be grossly negligent 
to perform this limit without explicitly knowing criteria, which the charge densities 
/ and g have to satisfy in order to guarantee that the limit /3 — is continuous and 
produces a finite result. 

This question can be analyzed with the help of Fourier transformation. If we 
use the symmetrical version of Fourier transformation according to (lA.ip and (lA.2p . 
the six-dimensional integrals (12. 4p and (12. 5p can be expressed as three-dimensional 
Fourier integrals via the following general expression introduced into electronic struc- 
ture calculations by Geller [581 Eqs. (1) and (2)]: 

r{r)g{r')h{r-r')drdr' = {27cf' j [f{p)Y-g{p)Hp)d'p. (2.6) 

If we choose h{r — r') = exp(— /5[r — r'])/\r — r'\, we only need the Fourier transform 
of the Yukawa potential (see for example [591 Eqs. (6.8) and (6.9)]), 

r exp(-ip.r) eM-M ^3^ ^ iVn)^ ^ ^ ^ ^ 
J r 
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to obtain: 



yif,9;P) = 4vr 1 [f{p)Y^^^g{p)d'p. (2.8) 

Fourier transformation maps L^(]R'^) onto L^(M^) in a one-to-one manner such that 
inner products are conserved pOl Theorem IX. 6 on p. 10]. Thus, u, f G L^(M^) 
imphes u,v & L^(]R^), and the coordinate and momentum space inner products (]B.3|) 
and (IB.gp . respectively, are identical: 



u{r)]*v{r)d^r = j [u{p)]* v{p) d?p (2.9) 



It follows at once from the Cauchy-Schwarz Inequality (12.21) that y{f,g', {3) is fi- 
nite if /(p)/[/3^ + e L2(M3)_ Unfortunately, this condition is 
not particularly helpful in practice since it is not so easy to find convenient coordi- 
nate representations for the functions /(p)/[/5^ + p^Y^"^ and g{p)/[l3'^ + p^Y^"^ (they 
have to be expressed as convolution integrals containing f{r) and g{r), respectively, 
multiplied by a modified Bessel function). 

Fortunately, a slightly more restrictive, but much more convenient condition on 
the charge densities / and g can be constructed easily. For all p e and for all 
/5 > 0, we have l/{(3'^+p'^) < l/f3\ Thus, 



1 

< — 



[fip)Y9ip) d'p 



(3>0. (2.10) 



By applying the Cauchy-Schwarz Inequality (12. 2p to the integral on the right-hand 
side we find that y{f,g',P) with P > is finite if f,g E L^(R^) or - since Fourier 
transformation is an isometric isomorphism of L^(R^) - if /, (7 G L^(M'^). 

If we perform the limit /3 — > in the Fourier integral (12. 8p . we formally obtain: 

C{f,9) = 47r I [fip)Y^g{p)d'p. (2.11) 

The momentum space integral on the right-hand side can be interpreted to be an 
inner product of the type of (IB. 11) that gives rise to a suitable Hilbert space. 

The Cauchy-Schwarz Inequality (12.21) now requires f{p)/p,g{p)/p G L^(R'^) to 
guarantee that C{f,g) is finite. Obviously, this is much more restrictive than the 
requirement f,ge L^jlR^) or even fip)/[P'^+p'^Y^'^,g{p)/[P'^+p'^Y^'^ E L'^{R^), which 
both guarantee that y{f,g]l3) is finite (alternative criteria, which also guarantee 
the existence of the Coulomb integrals C{f,g), are for instance formulated in [60} 
Example 3 (Sobolev's inequality) on p. 31] or [nU Section 4.3 (Hardy-Littlewood- 
Sobolev inequality)]). 

This example clearly shows that the limit /5 — in integrals involving the Yukawa 
potential exp(— /5r)/r is not necessarily continuous and does not always produce a 
finite result for arbitrary square integrable charge densities / and g. 

If we set /5 = in the Fourier transform (12.71) of the Yukawa potential, we formally 
obtain the Fourier transform of the Coulomb potential (see for example [621 Eq. (2) 
on p. 194]): 

(2.)-/^ / exp(-ip^^3^ ^ (2/vr)V^ 
J r p^ 
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There is a fundamental difference between the Fourier transforms (12 .7^ and (]2.12p . 
The Fourier integral in (12.71) is well defined and exists in the sense of classical analysis, 
whereas the Fourier integral in fl2.12l) diverges in the sense of classical analysis and 
becomes mathematically meaningful only if certain limiting or summation procedures 
are applied. Thus, the Fourier transform (12.71) of the Yukawa potential is a function 
in the ordinary sense. In contrast, the Fourier transform (12.121) of the Coulomb 
potential is a generalized function or distribution which is meaningful in suitably 
restricted functionals only. 

So, whenever we perform the limit — in integrals involving the Yukawa 
potential exp(— /?r)/r, we have to be cautious and take into account that this limit 
may be discontinuous and that it does not necessarily produce a finite result in the 
case of an essentially arbitrary integrand. It will become clear in later Sections 
that these problems are not restricted to integrals and that they can also occur in 
orthogonal expansions of the Yukawa potential. 

3 Divergent Orthogonal Expansions 

Hilbert spaces, whose basic features are reviewed in Appendix [Bl play a major role 
in various branches of mathematics and mathematical physics and in particular also 
in approximation theory. They also provide a rigorous mathematical framework for 
quantum mechanics. 

As discussed in Section H] or in more details in [HI Section 3], one-range addition 
theorems can also be viewed to be special approximation procedures: A function 
/(r ± r') belonging to suitable Hilbert space is expanded in terms of a complete and 
orthonormal function set in such a way that the two argument vectors r, r' G M? 
are separated. By construction, such an orthogonal expansion converges in the mean 
with respect to the norm of the corresponding Hilbert space. 

Therefore, it certainly makes sense to discuss the basic properties of orthogonal 
expansions in Hilbert spaces - including their power as well as their limitations - in 
a relatively detailed way. Hilbert spaces are linear vector spaces over the complex 
numbers equipped with an inner product (•!•) satisfying (IB.ip and a norm || ■ || sat- 
isfying (lB.2p . which has to be finite. A vector space with these properties is called a 
Hilbert space if it is complete with respect to its norm || ■ ||. 

Let us assume that / is an element of some Hilbert space 7i, and that the functions 
{(/9„}^Q are linearly independent and complete in 7i. Then, / can be approximated 
by finite linear combinations 

N 

fN = J^CfVn, NeNo. (3.1) 

n=0 

The coefficients Ci^^ are chosen in such a way that the mean square deviation 

ll/-/ivf = if-fNlf-fN) (3.2) 

becomes minimal. 
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The determination of the coefficients Cn in (13.11) by minimizing the mean square 
deviation (13. 2p only makes sense if both / and the functions {v?n}5^Lo normalizable 
according to ||/|| < oo and \\(Pn\\ < oo, respectively. Thus, / as well as the functions 
{(Pn}'i^=o have to belong to the Hilbert space TC. 

The finite approximation (13. ip converges to / as ^ cxd if the mean square 
deviation (13. 2p can be made as small as we like by increasing the summation limit 
N. In the case of convergence, it looks natural to assume that / possesses an infinite 
expansion 

oo 

/ = J2'^n^- (3.3) 

n=0 

in terms of the linearly independent and complete functions {(pn}'^=o with coefficients 

Cn = limAT^oo Ci!^\ 

Unfortunately, this is not true. In general, the coefficients Ci^^ in (13. ip do not 
only depend on n, /, and {(Pn}'^=o, but also on the summation limit N. It is not a 
priori clear whether the coefficients Cn^^ in (13. ip possess well defined limits Cn = 
limjv^oo C'n^^ or to put it differently, whether an infinite expansion of the type of 
(13.31) exists. Expansions of the type of (13. 3p may or may not exist. 

It is one of the central results of approximation theory that for functions / G 7i the 
mean square deviation (13. 2p becomes minimal if the functions {(pn}^=o are not only 
linearly independent and complete, but also orthonormal satisfying {fnl'Pn') = ^nn' 
for all indices n, n' G Nq, and if the coefficients are chosen according to Cn^^ = {^n\f) 
(see for example [IHl Theorem 9 on p. 51]). 

If the functions {(fn}'^=o ai's complete and orthonormal in Ti and if the expansion 
coefficients are chosen according to Ci^"* = (y9„|/), then the coefficients (v^nl/) in Jn 
do not depend on the truncation order A^. Thus, f E Ti possesses an infinite series 
expansion 

oo 

f = Y.{^n\f)^n (3.4) 

ri=0 

in terms of the complete and orthonormal function set {v?n}J^O' ^^"^ ^^^^ expansion 
converges in the mean with respect to the norm || ■ || of the Hilbert space Ti. 

This is all well known and described in countless books on functional analysis 
or approximation theory. In these books, it is always emphasized that orthogonal 
expansions of the type of (13. 4p are mathematically meaningful and converge in the 
mean if and only if / G 7i. 

However, in practical applications we are often confronted with functions that are 
not normalizable and thus do not belong to the corresponding Hilbert space Ti. In 
some cases it may be desirable to expand such a function / ^ 7i in terms of functions 
{<^ri}^0' t^^t complete and orthonormal in Ti. It is thus a practically relevant 
question whether and under which conditions the concept of orthogonal expansions 
in a Hilbert space Ti can be extended to functions f ^Ti. 

If this is indeed possible, we also have to analyze in which respect orthogonal 
expansions of a function / G 7i differ from those of a function f ^Ti. In particular. 
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we have to analyze whether and under which conditions orthogonal expansions / = 
Yl'^=oifn\f)fn can be used in inner products {f\g) with f and g ^ 7i. 

There is one obvious complication: li f & Ti., then it follows from the Cauchy- 
Schwarz Inequality (12.21) that the map g i— >■ {f\g) is bounded and thus continuous for 
all g E H. If, however, / ^ 7i, we have to take into account the Riesz Representation 
Theorem (see for example [171 Theorem 7.60 on p. 199]): 

For every continuous and thus bounded linear functional U: Ti ^ C there 
exists a unique u eH such that U{y) = {u\v) for all v ETi. 

Thus, the map g ^ {f\g) with f ^ Ti cannot be continuous and bounded for all 
g ETi. Consequently, for a given f ^H, there must be at least one g eH that yields 
an unbounded inner product {f\g)- 

Accordingly, the following discussion has to be limited to those g that yield 
for a given f ^ 7i bounded inner products {f\g). Therefore, we have to assume that 
g belongs to the subset JF c defined by 



For arbitrary f ^ H, it can happen that there is no 5^ G H satisfying K/lf?)] < 00, 
i.e., that JF is empty. In the following text, it will be assumed that this is not the 
case. However, it follows from the Riesz Representation Theorem that we cannot 
have T = TL, i.e., JF is either empty or a proper subset of Ti. 

A divergent orthogonal expansion / = Yl^={)i.Vn\f)Vn makes sense only if the 
inner products {^n\f) are finite for all finite indices n. If / G Ti, this is guaranteed 
by the Cauchy-Schwarz Inequality (12. 2p . but for / ^ 7i we have to assume explicitly 
that G J-" holds for all finite values of the index n. 

On the basis of these assumptions, it is at least formally possible to construct 
an orthogonal expansion / = Yl^=aiVn\f)Vn even if f ^ H. It is, however, not 
at all clear whether and in which sense the formal series expansion ^^o(v'"l/)v'" 
represents f ^ Ti. In particular, we have no a priori reason to assume that this 
expansion might converge in the mean according to the norm || ■ || of 7i, which would 
imply that Yln=o \ {^n\f)\^ < 00 holds. 

The divergence of the series Kv'nl/)!^ can be made plausible by analyzing 

the mean square deviation of the difference between a function / and its (possibly 
divergent) orthogonal expansion Yl'^=Qi.'^n\f)'{'n- 



T = {g\gen,f ^n,\if\g)\< 00}. 



(3.5) 



00 



00 00 



\\f-J2Mf)^^ 



'n 



) 



n=0 




00 



n=0 



00 



00 



00 



iVnlf)* iVmlf) iVmlVn) + ^ ('^'"1/) iflVm) 



m=0 n=0 m=0 



oo 




(3.6) 



m=0 
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If / G 7i, we have ||/||^ = (/|/) < oo, and it follows from ParsevaVs Equality (see 
for example [63], Eq. (II. 2) on p. 45]) 



that (13.61) vanishes. 

If, however, f ^ H, then = (/|/) diverges. In this case, the mean square 
deviation (13.61) can only vanish if J2m=o\i'^rn\f)\'^ also diverges. Of course, this is 
hand-waving and not a rigorous mathematical proof. In particular, it is not at all clear 
whether the mean square deviation (13.61) makes sense at all ii f ^ H. Nevertheless, 
this non-rigorous argument should suffice to convince even a skeptical reader that 
it would be overly optimistic to expect that the formal expansion '^'^=Q{'^n\f)'^n 
converges in the mean to something finite if / ^ H. 

We should also take into account that series expansions - or actually all approx- 
imation schemes - are practically useful only if they reproduce the essential features 
of the function they represent. At least in quantum mechanical bound state calcula- 
tions, the norm of a function is of considerable importance and should be preserved 
by an orthogonal expansion. Otherwise, a function / ^ 7i could be transformed to a 
function belonging to Ti simply be expanding / in terms of a complete and orthonor- 
mal function set. This would truly be a miraculous achievement with revolutionary 
and most likely highly undesirable consequences. 

Plausibihty arguments are no substitute for a rigorous proof. Such a proof can be 
formulated with the help of the Riesz-Fischer Theorem (see for example [171 Theorem 
7.43 on p. 191]): 

Let {ipn}'^=o ^ complete and orthonormal function set in a Hilbert space 
7i, and let cq, ci, C2, . . . he a sequence of numbers such that Yl'^=o 
converges. Then, the expansion X]^o ^nV^n converges in the mean to some 
w E H such that Cn = (ipnlw). 

This theorem shows that there is no sequence (ipolf), {ipi\f), {(p2\f), ... of inner 
products such that || J2'^=oi'Pn\f)fn\\'^ = J2n=o I(v5nl/)P < and f simultane- 
ously hold. Accordingly, f ^ H implies that J2'^=o('^n\f)fn diverges in the mean with 
respect to the norm of Ti. Thus, normalization is preserved by orthogonal expansions 
even ii f ^ H. 

This applies also to expansions in terms of Guseinov's functions ;j,\E'J^£(/3, r) defined 
by (ID.ip . As discussed in [HI p. 20], the Yukawa potential belongs for A; = 0, 1, 2, . . . 
to the weighted Hilbert space L^j.(M^) defined by fID.Sp . but not for k = —1 (see 
also [Ml P- 410]). Consequently, the expansion of the Yukawa potential in terms 
of Guseinov's functions converges in the mean for k = 0, 1,2, . . . and diverges for 
k = —1. In Appendix [EJ it is explicitly shown that this is indeed the case. 

We thus arrive at the conclusion that conventional Hilbert space theory and the 
related concept of convergent orthogonal expansions of the type of (13.41) in terms of 
complete and orthonormal functions {v^njJ^Lo only make sense if / G 7i. Divergent 



oo 




(3.7) 



n=0 
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expansions / = Yl'^=oi¥^n\f)'^n with f ^ TC are generalized functions in the sense of 
Schwartz [18] that can converge weakly when used in suitably restricted functionals. 

In the context of one-range addition theorems, which are to be used in multicenter 
integrals, it is essential that orthogonal expansions / = '^'^=Q{'^n\f)'^n can safely be 
used in inner products {f\g). If f,g E H, this can be shown by applying the Cauchy- 
Schwarz Inequality (12.21) to the inner product (/ — J2n=o('^n\f)'^n I d)- obtain: 



N N 



I (/ - Y.(^n\fWn I ^7) r < 11/ - J2(^n\fWnr M' ■ (3-8) 
n=0 n=0 

The assumption g E H implies WgW^ < oo, and the assumption f E Ti. implies that 
11/ — '^n=o('^n\f)fn\\'^ vauishes as oo. Accordingly, the right-hand side of (13. 8p 
vanishes as ^ oo. 

If g E TC but / ^ Ti, we cannot use the Cauchy-Schwarz Inequality (12.21) . More- 
over, the inner products {f\g) is not necessarily finite. As discussed above, the map 
g f— > {f\g) with f ^ Ti. cannot be continuous and bounded for all g E Ti. Thus, it is 
essential to assume that g belongs to the subset J-" C H defined by (13. 5p . 

Next, we have to formulate criteria, which guarantee that the divergent orthogonal 
expansion / = ^^o(v^"l/)v^" '^^^ safely be used in inner products {f\g) with g E J-'. 
Thus, we have to analyze under which conditions 

oo 

[f -Y.{Vn\f)^n\ g) = (3.9) 
n=0 

with f ^ Tt holds for suitable g E J-'. 

There is a very simple situation in which (13.91) is obviously valid. Let us assume 
that there is some g E that possesses a finite expansion in terms of the complete 
and orthonormal functions {v?n}5?Lo- 

M M 

g = ^TmV^m = (^|V5m)V5m, M E Nq . (3.10) 

m=0 m=0 

Since we always assume that the inner products {fn\f) are finite for all finite values 
of n, we obtain: 



{f -^i'^n\f)'^n\ g) = {f\g) - (^{¥^n\f)'^r^ 
n=0 n=0 

M oo M 

M M 
= ^ImiflVm) - ^{VmlfTlm = 0. (3.11) 



m=0 m=0 



Thus, for functions g E T satisfying (I3.10p . the divergent orthogonal expansion / 
^'^=Q{'^n\f)'-Pn produces the correct result in the inner product {f\g)- 
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As a mild generalization of (13.10P , let us now consider some g ^ T that possesses 
an infinite expansion in terms of the complete and orthonormal function {v?n}J^o- 

cc oo 
g = ^7m<y5m = ^{(Pm\g)^m- (3.12) 



m=0 m=0 



Since g E G Ti, this expansion converges in the mean. However, the convergence of 
this expansion alone does not suffice to guarantee that (13.91) is satisfied. The problem 
is that we are now confronted with infinite series that do not necessarily converge: 

oo oo 

(/ - '^{Vn\f)Vn \g) = (fig) - [^Mf)'^r. 
n=0 n=0 
oo oo 

= ^Imifl^m) - ^{^nlfY {Vn\g) 



m=0 n=0 
oo oo oo 



m=0 n=0 m=0 

oo oo 

7m{f\^m) - X^(</'n|/)*7n,- (3.13) 



m=0 n=0 



If I]m=o7m(/|v5m) and Zir=o(v'«l/)*7n both converge to {f\g), is satisfied, and 
the use of the divergent orthogonal expansion / = J2'^=oi'^n\f)'^n in the inner product 
{f\g) produces the correct result. 

The requirement, that the infinite series in (13.131) have to converge, makes it 
possible to characterize the subset JF c 7i defined by (13. 5p more precisely. The 
expansion g = Yl'^=oi'^n\g)'^n converges in the mean if the coefficients {y:>n\g) decay 
more rapidly than n"^/^ as n — > oo. Since / ^ 7i, the coefficients {^n\f) of the 
expansion / = X]^o('^nl/)'^" either decay less rapidly than n~-^/^ as n ^ oo or they 
may even diverge as n — oo. Thus, the coefficients {^Pn\g) in (I3.12p have to decay 
so fast as n ^ oo that the infinite series {f\g) = Yl'^=o(f\'^n){'^n\g) converges. This 
is certainly the case if the coefficients {f\ipn){'^n\g) decay more rapidly than 1/n as 
n ^ oo. 

Divergent orthogonal expansions / = J2'^=oi'^n\f)'^n with f ^ H possess the 
characteristic features of generalized functions in the sense of Schwartz [18] : Although 
divergent in the mean, such an expansion is meaningful in inner products {f\g) as 
long as g is restricted to the proper subset J-" C H defined by (13.51) . If g ^ J-", 
the divergent orthogonal expansion / = X]^^o(v^" l/)v^" cannot be used in the inner 
product {f\g) since the infinite series {f\g) = J2m=oif\'^rn){'^m\g) diverges. 

Conceptually, this situation very much resembles the theory of rigged Hilbert 
spaces or Gelfand triplets $ C C Here, is a Hilbert space, $ is a suitably 
restricted subset of H, and $^ is its dual space defined by the condition that inner 
product {u\v) with -u G $^ and f G $ remains finite. Loosely speaking, we may say 
that the more we restrict the subset $ C 7i, the larger its dual space becomes. 
A very readable account of rigged Hilbert spaces from the perspective of quantum 
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mechanics and their relationship with Dirac's bra and ket formalism can be found in 
the book by Ballentine [651 Chapter 1.4]. 

The insight, that divergent orthogonal expansions / = J2'i^^Q{fn\f)Vn with f 
are essentially generalized functions and can be used in a mathematical rigorous way 
in inner products {f\g) with g & J-', does not imply that all problems with the use 
of these divergent series are solved. The characterization of the subset JF c 7i is 
the crucial step that makes these divergent expansions mathematically meaningful 
in inner products {f\g)- But in real life applications, the characterization of JF may 
turn out be the most difficult problem that can occur in this context. 

4 One- Range Addition Theorems for the Coulomb 
Potential 

Let us assume that / belongs to the Hilbert space L^(M^) of square integrable func- 
tions defined by (IB.Sp and that the functions {'^ni{r)}n/,m are complete and or- 
thonormal in L^(]R^). As discussed in more details in [41 ^ Section 3], a one-range 
addition theorem for /(r ± r'), which converges in the mean with respect to the 
norm of L^(M^), can be constructed by expanding / in terms of the orthonormal 
functions {v25^£(T')}n,£,m: 

f{r ±r') = J2 CnAf: ±r') <,(r) , (4.1a) 

n£m 

CZif-^±r') = y"[<,(r)]*/(r±r')dV. (4.1b) 

The expansion fl4.ll) is indeed a one-range addition theorem, since the variables r 
and r' are completely separated: The dependence on r is entirely contained in the 
functions (/^^^(r), whereas r' occurs only in the expansion coefficients C^^{f; ±r') 
which are overlap integrals. 

If the overlap integrals Cl^g{f] ±r') can be expanded in terms of the functions 
(finei'^') according to 

Cn":.(/;±0 = E ^n4w(/;±X£'K), (4.2a) 

T-,\Z'{f;±) = I [<,,(r')]*C-(/;±r')dV, (4.2b) 
then the addition theorem fl4.ip assumes a completely symmetrical form: 

/(r±r') = ^n4w(/;±)<.W<,'(r'). (4.3) 

nim 
n'e'm' 

As is well known in approximation theory, a nontrivial weight function w{r) ^ 1 
can give more weight to those regions of space in which / is large, while deemphasizing 
the contribution from those regions in which / is small. Accordingly, the inclusion 
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of a suitable weight function t;;: M'^ ^ M+ can improve convergence. It is thus an 
obvious idea to construct one-range addition theorems that converge with respect to 
the norm of a weighted Hilbert space Iv^(M^) defined in (IB.Sp . If / G L^(M^), we 
can construct a one-range addition theorem by expanding /(r ±r') with respect to a 
function set {ipll!'£{r)}n/,m that is complete in L^(]R^) and orthonormal with respect 
to the modified inner product (IB.6P [HI Eqs. (3.5) - (3.7)))]. 

In the articles [3lll[5l[6l[7l[8l[9l[T0l[IIl[T2l[T3l[Tl[T5l^ 
[23l[2l[25l[26l[28l[29l[30l[ni[33[3l[^ Guseinov and coworkers derived and 

applied one-range addition theorems in connection with the special weight function 
w{r) = with k = —1, 0, 1, 2, . . . and used Guseinov's functions ^\E'™^(/9, r) defined 
by ( ID. Ill as expansion functions. This yields one-range addition theorems of the 
following general kind: 

fir ±r') = J2 kT^t'em'U: /3, ±) ^) XieW. r) , (4.4a) 

nlm 
n' I'm! 

.T::f,^'(/;/9,±) = / [.<.'(/9,^')]*(r')%C™,(/;A±r')dV, (4.4b) 

,C:^,(/;A±r') = 1 [,vl/™(/?,r)]V'=/(r±r')dV. (4.4c) 

If / G L^fc(R^), this addition theorem converges in the mean according to the norm 
dEl) of the weighted Hilbert space L^^i^). 

As discussed in more details in [41; Section 4], Guseinov derived in this way 
one-range addition theorems, for example for Slater-type functions with integral and 
nonintegral principal quantum numbers defined by (11. ip . For fixed fc = — 1,0,1,2,..., 
Guseinov's functions satisfy the orthogonality condition flD.2p and they are complete 



and orthonormal in the weighted Hilbert space L^j.(M^) defined by flD.5p . As long as 
the Slater-type functions, which are to be expanded, belong to the weighted Hilbert 
space L^s.(M^), this is a completely legitimate approach that leads to one-range addi- 
tion theorems for Slater-type functions which converge in the mean with respect to 
the norm (ID^ of L\(R3). 



The Coulomb potential plays a central role in electronic structure calculations, 
and the evaluation of inter-electronic repulsion integrals of the type of (12. 4p leads to 
formidable computational problems, in particular if the densities / and g in (12. 4p 
are two-center charge densities of the type of uir — A)v{r — -B) with A, S G M^. 
Therefore, it would be desirable to have a one-range addition theorem for the Coulomb 
potential. In [TTj, Guseinov tried to accomplish this by expanding l/|r — r'| in 
terms of functions fc\&^£(/?, r). Formally, Guseinov's approach leads to the following 
symmetrical one-range addition theorems: 

^ = E '^r;:f,7^'(/^)^^^^(/^'^)^<^'(/^'^'), (4.5a) 



r — r 



nlm 
n'£'m' 



.K'eM = J [kK'A(^,r')Yr'\C-/P,r')d'r', (4.5b) 
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,C^^,{py) = /[.^WO]*^:^dV. (4.5c) 

There are, however, some principal problem with these orthogonal expansions which 
Guseinov had either overlooked or ignored in his earlier work and which he still 
ignored in his most recent preprints [HI HSj |46], although I had emphasized their 
importance in [H]. In order to convince both Guseinov as well as other skeptical 
readers, I presented in Section [3] a detailed discussion of the properties of orthogonal 
expansions. The central features of these expansions can be summarized as follows: 

1. If / belongs to some Hilbert space 7i, then / can be expanded in terms of a 
function set {v^nj^o ^^^^ complete and orthonormal in ?i, and the expansion 
Yl^=Q{Vn\f)H^n converges to / in the mean with respect to the norm of Ti. 

2. If / ^ Ti, the formal orthogonal expansion / = X]^o(v^"l/)v^" diverges in the 
mean with respect to the norm of Ti. 

3. Nevertheless, such a divergent orthogonal expansion / = Yl'^=oi.Vn\f)Vn can 
produce meaningful results in inner products {f\g) as long as g is restricted to 
the proper subset JF c 7i defined by (13. 5p . 

The weighted Hilbert spaces L^t.(M^) with k = —1,0,1,2,... are based on the 
inner product flD.3P which involves an integration over the whole three-dimensional 



space with weight function w{r) = r^. The Coulomb potential l/\r — r'\ does not 
belong to any of the weighted Hilbert spaces L^fc(M^) which Guseinov implicitly used 
in his work. This implies that the one-range addition theorems (14.51) of the Coulomb 
potential diverge for k = —1, 0, 1, 2, . . . in the mean with respect to the norms flD.4p 
of the weighted Hilbert spaces L^s.(M^). 

Guseinov was not the first one who had derived a divergent expansion for the 
Coulomb potential in terms of a complete and orthonormal function set. Salmon, 
Birss, and Ruedenberg [66j derived a bipolar expansion of the Coulomb potential in 
terms of the Gaussian-type eigenfunctions of a three-dimensional isotropic harmonic 
oscillator which are complete and orthonormal in the Hilbert space L^(M^) of square 
integrable functions. However, Silverstone and Kay [67] demonstrated that this ex- 
pansion diverges, which was confirmed by Ruedenberg and Salmon [68] . Apparently, 
this observation was some kind of death sentence for the bipolar expansion of Salmon, 
Birss, and Ruedenberg [BB]. As far as I know, nobody has ever used this expansion. 

The divergence of the one-range addition theorems (14.51) for the Coulomb potential 
cannot be ignored. So, we are confronted with the question what we should do with 
these addition theorems. We could dismiss them as practically useless and ignore 
them, as it was done with the divergent bipolar expansion of Salmon, Birss, and 
Ruedenberg [66]. However, I think that this would be premature. Salmon, Birss, and 
Ruedenberg published their divergent bipolar expansion in 1968, and we now have 
a much better understanding of generalized functions and we also know much more 
about numerically efficient summation techniques which often can associate a finite 
value to a divergent series (see also the discussion in [HI pp. 24 - 25]). 
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As discussed in Section El divergent expansions of the type of (I4.5P can neverthe- 
less be practically useful as long as they are exclusively used in suitably restricted 
functionals. Thus, these functions are essentially generalized functions in the sense 
of Schwartz [Hj. Obviously, this offers new perspectives, but one must not forget 
that in relationships involving generalized functions one always has to be extremely 
careful about their domains of validity 

Consequently, it is not acceptable to proceed like Guseinov and to ignore the 
distributional nature of divergent one-range addition theorems of the type of (14. 5p 
and to treat them like ordinary orthogonal expansions that converge in the mean. 
It is absolutely essential to formulate regularity conditions for functionals - in our 
case multicenter integrals - and to take them into account. Otherwise, the use of 
distributional one-range addition theorems in multicenter integrals would be purely 
experimental. Virtually every outcome would be possible, depending on the other 
functions occurring in the integral. The fact that Guseinov apparently encountered 
no problems in his numerical examples does not prove anything. 

I do not expect that it will be easy to formulate the necessary regularity condi- 
tions. Firstly, the Coulomb potential is in spite of its apparent simplicity a relatively 
complicated mathematical object (see for example [611 Chapter 9]). Secondly, I fear 
that for every type of multicenter integral containing the Coulomb potential a new 
set of regularity conditions has to be formulated. 

In im p. 23]), a possible strategy based on orthogonal expansions of the charge 
densities / and g in Coulomb integrals C{f,g) defined by (12. 4p was sketched. It 
cannot be denied that this approach would be highly pedestrian, and more elegant 
and more powerful alternatives would be highly desirable. It seems that a lot of work 
remains to be done before distributional one-range addition theorems of the type of 
(14.51) can safely and effectively be applied in multicenter integrals. 

The idea of using distributional orthogonal expansions, which diverge in the mean 
and converge only weakly in suitably restricted functionals, is not new. In [69], I 
derived expansions of the plane wave exp(ip-r) in terms of complete orthonormal and 
biorthogonal function sets that converge only weakly. In some cases, these expansions 
simplify the evaluation of Fourier transforms, and they can also be used for the 
construction of one-range addition theorems (see [61] or [691 Section VII]). 

Expansions for the plane wave, that closely resemble those derived in [69], were 
also constructed by Guseinov [6| Eqs. (45) - (46)]. Guseinov, who did not mention 
[69] in [6], either overlooked or deliberately ignored the obvious fact that that the 
plane wave does not belong to any of the Hilbert spaces which he implicitly used. 
Accordingly, Guseinov's expansion diverge in the mean and can only converge weakly. 
Guseinov's oversight is hard to understand because he had cited [69] in several other 
articles [3lIIl^[T8l[l9lE7l[28l[32]. 

5 Guseinov's Rearranged Addition Theorems 

As discussed in Section [3l orthogonal expansions play a central role in Hilbert spaces 
and also in approximation theory. In contrast, nonorthogonal expansions are largely 
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ignored. Of course, there are many good reasons for this neglect. In the context 
of one-range addition theorems, the most important evidence speaking against the 
use of nonorthogonal function is the following well established fact: If a function 
set {9?n}^o is only complete in a given Hilbert space 7i, but not orthogonal, then 
it is general only possible to construct finite approximations to / G 7i of the type 
of (13.11) by minimizing the mean square deviation (13.21) . but the existence of formal 
expansions of the type of (13.31) in terms of nonorthogonal functions is not guaranteed: 
Thus, nonorthogonal expansions may or may not exist. This fact is documented quite 
extensively in the mathematical literature (see for example |^ Theorem 10 on p. 54] 
or [501 Section 1.4]) or also in the literature on electronic structure calculations [511 
[521 [531 [5^ [551 [56]). Horrifying examples of pathologies of nonorthogonal expansions 
can be found in [SSI Section III. I]. 

Of course, there are situations in which nonorthogonal expansions offer compu- 
tational advantages (see fore example the discussion in t2D|)- However, in the vast 
majority of all cases, orthogonal expansions have clearly superior properties. Conse- 
quently, one should not voluntarily abandon the highly useful feature of orthogonality 
unless there are truly compelling reasons. 

In [Tl[3], Guseinov derived one-range addition theorems for Slater-type functions 
with integral and nonintegral principal quantum numbers of the type of (14. ip by 
expanding them in terms of his complete and orthonormal functions ;i,\E'™^(/5, r). 
As long as the Slater-type functions belong to the weighted Hilbert spaces L^fc(M^) 
defined by flD.5l) with k = —1,0,1,2,..., these addition theorem converge in the 
mean with respect to the norms flD.4|) of these Hilbert spaces. 

For reasons, which I do not really understand, Guseinov considered it to be ad- 
vantageous to replace in his one-range addition theorems for Slater-type functions his 
complete and orthonormal functions ^\E'^£(/5, r) by nonorthogonal Slater-type func- 
tions with integral principal quantum numbers via [HI Eq. (6.4)] 

and to rearrange the order of summations of the resulting expansions. Guseinov 
constructed in this way expansions of Slater-type functions Xn^lH^^^ ^ with in 
general nonintegral principal quantum numbers G R \ N in terms of Slater-type 
functions x^ii.f^i'"^) with integral principal n G N quantum numbers located at a 
different center (see also [H] Section 6]). 

As is well known. Slater-type functions are complete in all Hilbert space implicitly 
used by Guseinov (for a proof, see [HH Section 4]), but not orthogonal. In view of 
the principal problems mentioned above, it is therefore not at all clear whether Gu- 
seinov's rearranged addition theorems are mathematically meaningful. Accordingly, 
I claimed in [HI Section 6]) that Guseinov's rearrangements of his one-range addi- 
tion theorems, which are expansions in terms of his functions ^\&^£(/3, r) and thus 
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ultimately expansions in terms of generalized Laguerre polynomials L^^^^_^^^^(2/5r), 
are dangerous and potentially disastrous and that their validity has to be checked. 

In [m p. 7], Guseinov disagreed and claimed this his Eq. (3.11) - a finite nested 
sum containing his functions f^^if^^^P, r) on the left-hand side and Slater-type func- 
tions X^iiP, r) with integral principal quantum numbers on the right-hand side - 
proves the validity and mathematical soundness of his approach. 

I do not question the validity of Guseinov's Eq. (3.11), but I very much disagree 
with Guseinov's conclusion that his Eq. (3.11) proves the validity of his rearrange- 
ments. The problem with Guseinov's reasoning is that he does not distinguish care- 
fully between rearrangements of finite and infinite sums. Obviously, a finite sum 



TV 



F^{x) = 5^ A(r)L(f)(x), NeNo, (5.2) 



n=0 



of generalized Laguerre polynomials lI^\x) multiplied by purely numerical coeffi- 



cients aI"^ can always be rearranged. If we insert the explicit expression (IC.ip of the 
generalized Laguerre polynomials into (15. 2p and rearrange the order of summations, 
we obtain 

F„W = f |;(^±i:±i).,(t. (5.3) 

But if we now perform in the finite sum (15. 2p the limit N ^ oo and consider instead 
the rearrangement of the infinite series 

oo 

F{x) = hm Fr,{x) = Ai") 4") (x) , (5.4) 

n=0 

the situation is much more complicated and many things can go wrong. Formally, a 
rearrangement of F{x) yields the following power series in x: 

^w-E^E^^^^^'-t. (5.5) 

This power series for F{x) makes sense if and only if the inner series on the right-hand 
side of (15. 5p converges for every z/ G Nq. In addition, the inner series in has to 
produce values that do not increase too strongly with increasing z/, because otherwise 
the right-hand side of (15. 5p diverges for every |x| > 0. In the case of an essentially 
arbitrary function F[x), these two conditions are not necessarily satisfied. It is also 
easy to show that the convergence of the Laguerre expansion for F{x) with respect 
to the norm of the Laguerre-type Hilbert space L^_^^„(M_|_) defined by (10.30 does not 
imply the convergence of the inner series over on the right-hand side of (15.50 for 
every z/ G Nq. 

Special attention deserves the case that F(x) is not an analytic function in the 
sense of complex analysis at the expansion point x = 0. In this case, the rearranged 
power series on right-hand side of (15. 5p cannot exist because all but a finite number 
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of series coefficients are infinite, even if the Laguerre expansion for F{x) exists and 
converges in the mean with respect to the norm of Lg_^^cj(M+). 

Accordingly, I see no reason to aher my assessment ^Ti p. 18] that Guseinov's 
rearrangements are dangerous and potentially disastrous, and that their validity must 
be checked explicitly. 

The rearrangements of the finite sum F]\f{z) or of the infinite series F{z) yielding 
( 15.31) and (15. 5p . respectively, are special cases of the rearrangements of double series 
12m=o X]^o ^rn,n- This is an old and extensively studied topic in the theory of infinite 
series. The most detailed treatment, which I am aware of, can be found in the book 
by Bromwich ^TJj Chapter V]. Loosely speaking, the rearrangement of such a double 
series is safe if the double series converges absolutely. In the case of expansions in 
terms of orthogonal polynomials, we cannot tacitly assume absolute converge. 

So, we have good reason to assume that Guseinov's rearrangements of one-range 
addition theorems are dangerous and potentially disastrous. Of course, it is not 
satisfactory if we only know that a given mathematical manipulation is dangerous. 
Instead, we would like to know with certainty whether this operation is legitimate or 
not. 

Unfortunately, one-range additions theorems for exponentially decaying functions 
are fairly complicated mathematical objects, and explicit proofs of their convergence 
and/or divergence are very difficult. Most likely, such an investigation would be a 
nontrivial research problem in its own right. Since I am convinced that Guseinov's 
rearrangements are basically a bad idea, I saw no point in spending too much time 
and effort. Therefore, I looked for simpler alternatives to a detailed convergence 
analysis, even if these alternatives would not answer all questions of interest. 

As shown in [1T| Section 6], valuable insight can in some cases be gained by con- 
sidering not the complicated one-range addition theorems themselves, but their much 
simpler one-center limits. Let us therefore assume that we succeeded in constructing 
a one-range addition theorem of the type of (14.41) for some function f{r± r') by 
expanding it in terms of Guseinov's functions. If we now consider its one-center limit 
by setting r' = 0, our addition theorem must simplify to yield an expansion of /(r) 
in terms of Guseinov's functions: 

fir) = ,^„":,(/5; /) ,KAf^^ ^) ' (5-6a) 

n£m 

k^uMJ) = I [fcVl/™(/?,r)]*r^/(r)dr. (5.6b) 

Under fortunate circumstances, the mathematical nature of such an identity allows 
conclusions about the legitimacy of Guseinov's rearrangements. 

Let us now assume that we succeeded in deriving a one-range addition theorem 
of the type of (14.41) by expanding Slater-type functions Xn liPy r±r') with in general 
nonintegral principal quantum numbers G M \ N in terms of Guseinov's functions 
j,\E'^£(/3, r) with equal scaling parameters /3 > (see for example [HI Eq. (6.1)] with 
f3 = 'j). If we now set r' = in this addition theorem, it must simplify to yield the 
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following expansion of liP^^) terms of Guseinov's functions: 

X^Af3,r) = T{N + L + k + 2) 



oo 



2^ y, 1 oM .ll/2 k'^y+L+lMI^^^) 



A^gR\N, /3>0, = -1,0,1,2,... . (5.7) 



If G N and > L + 1, the infinite series on the right-hand side terminates because 
of the Pochhammer symbol (— + L + 

It is easy to show that the expansion (15.71) in terms of Guseinov's function is a 
special case of the following expansion |^ Eq. (5.17)] which expresses a nonintegral 
power with yU G M \ Nq as an infinite series of generalized Laguerre polynomials: 

^„ _ r(M:a + i) ^ 

r(a + 1) (a + l)n 

/xGR\No, Re(/i + a)>-l, Re(a)>-1. (5.8) 



If we set jj = m with m G No, the infinite series on the right-hand side terminates 
because of the Pochhammer symbol (— m)„. 

In [4T1 Eqs. (6.9) - (6.11)], it was shown that it is not possible to transform the 
Laguerre expansion (15.81) for x^ with G M \ No to a power series in x by inserting 
the explicit expression ( IC.ll) of the generalized Laguerre polynomials. Interchanging 
the order of the nested summations yields a formal power series in x Eq. (6.9)]. 
Superficially, this looks like success. However, the coefficients of this power series can 
be expressed as hypergeometric series iFq which are for all but a finite number of 
indices infinite PT, Eq. (6.11)]. 

Of course, this failure is not really surprising: The general power function z'^ with 
z G C and /i G C \ No is not analytic at z = in the sense of complex analysis. For 
H = m with m G Nq, Taylor expansion of z"^ around z = is, however, legitimate 
and yields the trivial identity z"^ = z"^. 

Thus, we can conclude that in the case of equal scaling parameters /3 > the 
one-center limit r' = of Guseinov's rearranged addition theorem for Slater-type 
functions Xjvl(/^5 does not exist if the principal quantum number is nonin- 

tegral, A^ G M \ N. 

Let us now assume that we succeeded in deriving a one-range addition theorem 
of the type of (14. 4p by expanding Slater-type functions Xn l^P^ r ±r') with in general 
nonintegral principal quantum numbers A^ G M \ N in terms of Guseinov's functions 
^,^^^£(7, r) with different scaling parameters /3 7^ 7 > (see for example [HI Eq. 
(6.1)]). If we now set r' = in this addition theorem, it must simplify to yield 
the following expansion of xjlf^(/5, r) in terms of Guseinov's functions with different 
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scaling parameter 7 7^ /3 > 0: 



X 

iy=0 



{u + 2L + k + 2)\ ^ 
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X -z^, + + A; + 2; 2L + + 3; 



/5 + 7. 

A^gM\N, /3,7>0. (5.9) 

If we set 7 = /3, we of course obtain (15.71) . 

It is easy to show that (15.91) is a special case of the following expansion [HI Eq. 
(6.12)]: 



r(a + /i + l) 
r(a + l) 

n=0 ^ 

/ieM\No, Re(/i + a)>-l, mg(-oo,1/2). (5.10) 



The condition —00 < -u < 1/2 is necessary to guarantee that this expansion converges 
in the mean with respect to the norm of the weighted Hilbert space L^.^j^^ (IR+). For 
M = 0, (15.101) simplifies to give (15. 8p . 

If we insert the explicit expression fIC.ip of the generalized Laguerre polynomials 



into (15.101) and interchange the order of summations, we also obtain a formal power 
series in x. Unfortunately, an analysis of the resulting power series becomes very 
difficult because of the terminating Gaussian hypergeometric series 2F1 in (I5.10p 
(probably, an analysis of the behavior of this 2-F1 as n — > 00 would be a nontrivial 
research project in its own right). However, we can argue that the function exp('U2;) 
with /i, M, z G C is only analytic at 2 = in the sense of complex analysis if is a 
nonnegative integer, fi = m with m G Nq, yielding exp{uz) = 'u"-2™'*'"/n!. If 

fj, is nonintegral, yU G C \ No, a power series expansion of of exp(uz) around z = 
does not exist. 

Thus, also for different scaling parameters /3 7^ 7, the one-center limit r' = of 
the rearranged addition theorems for x]!^^(/5, r ± r') does not exist if the principal 
quantum number is nonintegral, G M \ N. 

Apparently, Guseinov deliberately ignores even now the fact that the Laguerre 
expansions (15. 8p for x'^ and (I5.10p for x'^e"^ cannot be transformed to power series 
expansions in x if /i G M \ No, although this had been emphasized in pp. 18 - 
19]. In [ini Eqs. (5) - (6)], Guseinov expanded Slater-type functions Xml^P^^^) with 
nonintegral principal quantum numbers as an infinite series of Slater-type functions 
X^^(/5, r) with integral principal quantum numbers, although these expansions do not 
exist since their terms are for all but a finite number of indices infinite. In ^5l Eqs. 
(7) - (10)], Guseinov tried to resell essentially the same nonexisting expansion for 
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what he calls Coulomb- Yukawa like correlated interaction potentials, which are apart 
from a different normalization nothing but special Slater-type functions. 

From a mathematical point of view, a one-range addition theorem for a function 
/(r ± r') is a mapping x — C. In my opinion, one-range addition theo- 
rems have the highly advantageous feature that they provide a unique infinite series 
representation of f{r± r') with separated variables r and r' that is valid for the 
whole argument set M.^ x M^. If we accept this premise, then we have to conclude 
that Guseinov's manipulations, which produced his rearranged addition theorems for 
Slater- type functions Xn r ±r'), are at least in the case of nonintegral principal 
quantum numbers G M \ N a complete failure. 

The analysis of the one-center limits of rearranged one-range addition theorems 
provides valuable insight in the case of Slater-type functions with nonintegral prin- 
cipal quantum numbers, but it does not answer all questions. In particular, my 
nonanalyticity argument allows no conclusions about the validity of Guseinov's re- 
arrangements in the case of Slater-type functions with integral principal quantum 
numbers. Another interesting but open question is whether Guseinov's rearrange- 
ments are in the case nonintegral principal quantum numbers invalid for the whole 
argument set x or whether they are invalid only in the one-center limit. This 
is a practically very relevant question. If only the one-center limit is invalid, then it 
would be conceivable that Guseinov's rearranged one-range addition theorems might 
be mathematically meaningful or possibly even numerically useful in a restricted sense 
as approximations, although they do not exist for the whole argument set x M^. 
This has to be investigated. 

These examples show that the situation is much more complicated than originally 
anticipated by Guseinov. Obviously, a lot of work remains to be done before we 
can claim with some confidence that we understand the subtleties of Guseinov's 
rearrangements sufficiently well. It should also be clear that the burden of proof lies 
in all cases with Guseinov. 

Nevertheless, I do not think that it would be a good idea to invest too much time 
and effort into an analysis of these most likely very difficult open questions. In my 
opinion, it is simply a bad idea to construct one-range addition theorems that use 
nonorthogonal functions as expansion functions. It would be much better to focus 
on those one-range addition theorems that are expansions in terms of complete and 
orthonormal function sets. 

The principal superiority of orthogonal expansion functions becomes particularly 
evident in the case of one-range addition theorems of the type of f l4.5p that do not 
converge in the mean with respect to the norm of an appropriate Hilbert space, but 
only weakly in the sense of generalized functions in suitably restricted functionals. 

For example, in Sectional analyzed under which conditions inner products {f\g) 
with f ^ Ti and g & Ti are mathematically meaningful and whether the divergent 
orthogonal expansion / = Yl^=oi.'^n\f)'{'n nan be used in these inner products. I 
showed that if g belongs to the subset JF c 7i defined by (13.51) . then the expansion 
{f\9) = Y.n=oif\v^n){^n\g) couverges if the expansion coefficients {f\ipn){(Pn\g) decay 
more rapidly than l/ri as n — ^ cxo. 

Ignoring all questions of convergence or existence, let us now assume that both / 
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and g can be expanded at least formally in terms of a complete, but nonorthogonal 
function set {ipn}'^=o- 

oo 

/ = J^F^V^^, (5.11) 

m=0 

oo 

g = Y.G^fiJn. (5.12) 

n=0 

If we insert these expansions into the inner product {f\g), we formally obtain the 
following double series: 

oo oo 

= E E [^^TG^f i^rnl^n) • (5.13) 
m=0 n=0 

There can be no doubt that it would be much harder to formulate convergence criteria 
for this complicated double series than for the comparatively simple series {f\g) = 
'^m=oif\^rn){^m\g) which wc obtaiu if we expand / and g in terms of a complete 
and orthonormal function set {v^njJ^Lo- 

This simple example shows that it is a highly dubious idea to expand generalized 
functions in the sense of Schwartz [58] in terms of nonorthogonal function sets. At 
best, we would be confronted with nontrivial technical problems. 



6 Summary and Conclusions 

To some extent, divergent series are the dominant theme of this Reply. I am fully 
aware that divergent series have been and to some extend still are a highly contro- 
versial topic. There are still many who thoroughly dislike divergent series and think 
that they should be banned from the realm of rigorous mathematics. In their opin- 
ion, divergent series should at best be considered to be some kind of mathematical 
pornography. In addition, there are many others who - either because of ignorance 
or because of over- confidence - wrongly believe that divergent series cannot occur in 
their work and who thus tend to ignore questions of convergence. 

It is now widely accepted that divergent series play a very useful role. They 
are indispensable tools in mathematics and in particular also in the mathematical 
treatment of scientific problems. Skeptical readers, who still prefer to ignore divergent 
series, should search Google Scholar (http://scholar.google.coni/) for "divergent 
series" or for related topics. They will be surprised by the large number of applications 
of divergent series in different scientific disciplines. 

There are principal differences between the divergent orthogonal expansions con- 
sidered in this Reply and the more familiar divergent power series, which for instance 
occur abundantly in quantum mechanical perturbation expansions or as asymptotic 
expansions for special functions. Divergent power series can be used for the numerical 
evaluation of the function they represent: With the help of suitable summation tech- 
niques as for instance Borel summation. Fade approximants, or nonlinear sequence 
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transformations it is frequently possible to associate a finite value to a divergent 
power series. 

In contrast, it is not intended to use the divergent orthogonal expansions of this 
Reply for the direct numerical evaluation of the function they represent. We only 
want to use these expansions in suitable functionals - typically multicenter inte- 
grals - because we hope for some formal simplifications. Actually, this applies to all 
one-range addition theorems: They are only intermediate results which ultimately 
produce series expansions for multicenter integrals. 

The use of convergent expansions in integrals has the undeniable advantage that 
normally only comparatively mild assumptions are needed to guarantee that integra- 
tion and summation can be interchanged and that the resulting expansions converge. 
Nevertheless, the use of convergent expansions in integrals is to some extent a luxury 
and not strictly necessary. We are free to use a divergent expansion in an integral 
and interchange integration and summation if we can guarantee that the resulting 
expansion converges to the correct result. 

Obviously, such an approach gives us additional possibilities, but it would be 
naive to expect a free lunch: It is grossly negligent to use divergent series in in- 
tegrals without explicitly knowing criteria of manageable complexity that guarantee 
the convergence of the resulting expansions. This is probably the most serious flaw of 
Guseinov's work on one-range addition theorems. Since he is apparently completely 
unaware of the fact that divergent orthogonal expansions occur in his work, he has no 
reason to think about additional criteria which could justify the use of his divergent 
expansions in multicenter integrals. 

In Abstract], Guseinov claims that all his formulas were numerically tested, 
but this does not prove anything. A much more profound understanding of these 
divergent expansions and their domains of validity is needed, before they could be 
applied safely and in a mathematically rigorous way. Otherwise, the use of Guseinov's 
divergent one-range addition theorems in multicenter integrals would be purely ex- 
perimental. 

Misconceptions about divergent series are also the core of Guseinov's problems 
with his rearrangements of one-range addition theorems discussed in Section |5l It 
seems that Guseinov wrongly believes that it is always safe to rearrange the order of 
summations of double series ^^=0 X]^o 

In recent years, Guseinov and coworkers were able to publish a remarkably large 
number of articles on one-range addition theorems [Il[2l[3llll[5l[6l[7l[8l[9l [101 ttH [121 

[I31[H[ia[ISl[ni[ianHlS[2IlE2lE31[2aE3[2Sl[2ZlEHl[2^ 

[571 [551 [3S1 [ID]- This fact and the dubious quality of these articles raise obvious 
doubts about the quality of our scientific publication system based on anonymous 
peer-refereeing, and also on the competence of Guseinov's referees who also failed to 
understand the mathematical subtleties of one-range addition theorems. 

Of course, Guseinov disagreed with this conclusion which was first expressed in 
HB p. 5]. In [H p. 3], he stated: 

The respectable referees very well understand and examined the published 
by Guseinov and his coworkers in the years from 1978 to 2006 papers on 
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one-range addition theorems. 

Unfortunately, I am not so optimistic. But it would be unfair to blame exclusively 
Guseinov's referees. Refereeing Guseinov's manuscripts is certainly not easy. It is 
Guseinov's trademark to produces a large number of short and largely overlapping 
articles on essentially the same topic. This makes it very hard even for a very com- 
petent referee not to get lost in Guseinov's flood of publications and to keep track 
of Guseinov's truly new results. Moreover, as I know from my own experience as 
a referee, there is always the temptation to be less critical in the case of a short 
manuscript than in the case of a (very) long manuscript. 

In my opinion, part of the problem are short articles. While there can be no doubt 
that short articles are well suited to present new experimental or computational 
results, they are basically unsuited for predominantly theoretical or mathematical 
topics. 

For example, in a theoretical article on multicenter integrals, it is first necessary 
to provide a usually (very) long list of special functions and other abbreviations 
and conventions. Then, it is necessary to give a compact, but hopefully sufficiently 
comprehensive description of the mathematical techniques, which are to be employed. 
To do these things in a reasonable and for the reader beneficial way, we need at least 
a few pages before we can start with the derivation and description of new results. 

If we nevertheless insist on writing (very) short articles, we can either shrink 
the in my opinion very important introductory part to an absolute minimum, or we 
can try to split the new results into numerous small pieces. Either alternative is 
undesirable: If we choose the first alternative, essential background information may 
be lacking and readability will most likely suffer quite a bit, and if we choose the 
latter alternative, we have to write numerous articles on essentially the same topic 
that contain virtually nothing new. 

Of course, a compromise would also be possible: One could write a large number 
of articles with highly condensed and thus more or less incomprehensible introductory 
parts, that also present at best infinitesimal increments of insight. 

The problems with short articles can be demonstrated convincingly by Guseinov's 
recent reprint [46J. Guseinov's only new result, which I could detect, are his Eqs. 
(6) - (7), which express a Slater-type function xjif^(/?,r) with nonintegral principal 
quantum numbers as an infinite series of his functions fc^E'™^(7, r) and which corre- 
sponds to (15. 9p . In his Eqs. (8) - (9), Guseinov tried to resell his Eqs. (6) - (7) as 
an expansion for what he calls Coulomb- Yukawa like correlated interaction poten- 
tials, which are apart from a different normalization nothing but special Slater-type 
functions. 

It is fairly easy to derive Guseinov's new expansion. We only have to combine some 
well known properties of generalized Laguerre polynomials with a formula from the 
book by Gradshteyn and Rhyzhik [721 Eq. (7.414.7) on p. 850]. Of course, Guseinov's 
new expansion can be published, but one may wonder whether Guseinov's new series 
expansion alone justifies a new article. A change of the editorial policy.of scientific 
journals with respect to short articles on predominantly mathematical topics might 
be helpful. 
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A Terminology and Definitions 

For the set of positive integers, I write N = {1,2,3,...}, and for the set of non- 
negative integers, I write No = {0, 1,2,.. .}. The real and complex numbers and the 
set of three-dimensional vectors with real components are denoted by M, C, and M^, 
respectively. M+ is the set of real numbers > 0. 

Fourier transformation is used in its symmetrical form, i.e., a function f : M.^ ^ C 
and its Fourier transform / are connected by the integrals 

f{p) = i2n)-'/'l e-P-V(r)dV, (A.l) 
/(r) = (27r)-3/2| e^-P/(p)d3p, (A.2) 

B Hilbert Spaces 

Let V be a vector space over the complex numbers C that possesses an inner product 
(•|-) : V X V — *• C, satisfying for all u,v,w eV and for all a G C [631 P- 36] 

{u\u) > 0, (B.la) 

{u\u) = u = 0, (B.lb) 

{u\v + w) = {u\v) + {u\w) , (B.lc) 

{u\av) = a{u\v) , (B.ld) 

{u\v) = {v\u)* . (B.le) 

Another essential concept is the norm || ■ || : V ^ M of the elements of a vector 
space V, satisfying for all u,v E V and for all a G C [63, p. 8] 

||n|| > 0, (B.2a) 

= ^ u = 0, (B.2b) 

\\<yu\\ = \a\ \\u\\ , (B.2c) 

\\u + f II < ||n|| + ||i;|| . (B.2d) 

Obviously, ||m|| = ^/{u\u) with m G V is a norm satisfying these conditions. 

A vector space V over the complex numbers C is called a Hilbert space, if it 
possesses an inner product (-I-) satisfying (IB. II) . and if V is complete with respect 
to the norm defined by ||u|| = ■>/ {u\u). Completeness implies that every Cauchy 
sequence in V converges with respect to this norm to an element of V. 

In bound-state electronic structure calculations, we have to take into account 
Born's statistical interpretation of the wave function. Thus, an obvious inner prod- 
uct for functions f,g:M.^^C, that can be used as basis functions in atomic and 
molecular electronic structure calculations, can be defined according to 

(119)2 = J [f{r)Yg{r)d'r. (B.3) 
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As usual, the integration extends over the whole M^. 

On the basis of the inner product (IB.Sp . the norm of a function /: — * C is 
defined according to 



= (B.4) 

The Hubert space L^(M'^) of square integrable functions is defined via the norm 
( IB.4p according to 

L2(m3) = |/:M3^c| J |/(r)|2d-V < cx)} 

= {/:M3^C|||/||2<oo}. (B.5) 

The formalism of Hilbert spaces can be generalized to include weight functions. If 
w : M'^ ^ is a suitable positive weight function, we define the inner product with 
respect to the weight function w for functions /, (7 : M"' — C according to 

{f\9U2 = I [f{r)Yw{r)g{r)d'r. (B.6) 

As in (IB.Sp . the integration extends over the whole M^. It is easy to show that {f\g)w,2 
with w{r) > is indeed an inner product satisfyng (IB.ip . 

On the basis of the inner product (lB.6p . the norm of a function /: M'^ — C with 
respect to the weight function w is defined according to 



w,2 



(/l/W- (B.7) 



The Hilbert space L'^{M.^) of square integrable functions with respect to the weight 
function w is defined via the norm (IB.7P according to 



L2(r3) = |/:R3^C J w(r) |/(r)|2dV < 00^ 

= {/:R3^C|||/|U,2<oo}. (B.8) 

It is not necessary to use the coordinate representation for the definition of the 
Hilbert spaces L^(M3^. Instead, the momentum representation can also be used. This 
is a consequence of the well-known fact that Fourier transformation defined via (lA.ip 
and (lA.2p maps L^(R3) onto L^(R3) in a one-to-one manner such that inner products 
are conserved [SUl Theorem IX.6 on p. 10]. Thus, f,g & L^(R3) implies that the 
Fourier transforms f{p) and g{p) are also elements of L^(R^). In addition, the inner 
product fIB.Sp satisfies 



if\9)2 = I [f{p)Y9{p)d'p. (B.9) 
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C Laguerre polynomials 

The generalized Laguerre polynomials L^n\x) possess the following explicit expres- 
sions [731 P- 240]: 

^"^^''"iFi(-n;a + l;a;). (C.lb) 



n! 



The generalized Laguerre polynomials satisfy for Re(Q;) > —1 and m, n G Nq the 
following orthogonality relationship [73], p. 241]: 

r x-^-^Lt\x)L^:\x)d.x = Ii^±^^5_. (C.2) 
Jo n\ 

The completeness of the generalized Laguerre polynomials in the weighted Hilbert 
space 



|/:M+^C ^ e-^x"|/(x)|2da; < oo} (C.3) 



is a classic result of mathematical analysis (see for example [501 P- 33], [74, pp. 349 
351], or [751 PP- 235 - 238]). 



D Guseinov's Function Set 

In [21 Eq. (1)], Guseinov introduced a fairly large class of complete and orthonormal 
functions which can be expressed as follows: 



(2/3)^+=^(n-£ - 1)! ^ 
{n + £ + k + l)\ 



^.pr 2.(f^+_Y'^(2/5r) yri^Pr) . (D.l) 



The indices satisfy n e N, A; = -1, 0, 1, 2, . . . , £ G No < - 1, < m < £, and the 
scaling parameter [3 is positive. 

Guseinov's functions satisfy the orthonormality relationship (compare also [21 Eq. 
(4)]) 

j [,*-,(/?, r)]V%vi/-',,(/?,r)dV = 5nn'5u'5^m'. (D.2) 

Accordingly, Guseinov's functions are a complete and orthonormal set in the weighted 
Hilbert space L^fe(M^) with k = —1, 0, 1, 2, ... , which is defined via the inner product 

if\9U,2 = j[f{r)Yr^g{r)d'r (D.3) 



32 



Ernst Joachim Weniger 



and the norm 

\\f\\r^2 = \J{f\9)rK2 (D.4) 

according to 

Ll,{M?) = |/:M=^^C j r'=|/(r)|2dV < 00} 

= {/:M3^C|||/||,.2<oo}. (D.5) 

E Expansion of the Yukawa Potential 

It is relatively easy to construct an expansion of the Yukawa potential in terms 
of Guseinov's function j;.\E'™£(/9, r). The Yukawa potential is a special Slater-type 
function satisfying 

(47r)^/Vx%(Ar). (E.l) 



r 

Setting Ar = L = M = Oin ([52D yields: 



[27r/(2/3)^'+i] T{k + 2)Y, .^°+i,o(/?> r) . (E.2) 

This expansion converges in the mean with respect to the norm (1D.4P of the weighted 
Hilbert space L?i,{E?) with k = —1, 0, 1, 2, . . . if the squares of the coefficients on the 
right-hand side decay more rapidly than 1/z/ as — 00. 

The behavior of the coefficients in (lE.2p can be analyzed with the help of the 
following asymptotic expression for the ratio of two gamma functions flEl Eq. (6.1.47) 
on p. 257]: 

= + O , .^00. (E.3) 

We then obtain the following asymptotic estimates for the zz-dependent part of the 
square of the coefficients on the right-hand side of (]E.2p : 

2 

^' u-^-^ + 0{u-''-^) , z/->oo.(E.4) 



[{u + k + 2)\u^] 



1/2 



[u + k + 2) 



Thus, the expansion (]E.2|) converges in the mean with respect to the norm (1D.4|) 
of the weighted Hilbert space L\(M^) for k = 0,1,2... and diverges for k = — 1. 



This is in agreement with the fact that the Yukawa potential belongs to the weighted 
Hilbert space L^fe(M^) for = 0, 1, 2 . . . , but not for k = —1. 
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